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Abstract. We study the time evolution of quantum one-dimensional gapless systems evolving from 
initial states with a domain-wall. We generalize the path-integral imaginary time approach that 
together with boundary conformal field theory allows to derive the time and space dependence of 
general correlation functions. The latter are explicitly obtained for the Ising universality class, and the 
typical behavior of one- and two-point functions is derived for the general case. Possible connections 
with the stochastic Loewner evolution are discussed and explicit results for one-point time dependent 
averages are obtained for generic k for boundary conditions corresponding to SLE. We use this set of 
results to predict the time evolution of the entanglement entropy and obtain the universal constant 
shift due to the presence of a domain wall in the initial state. 
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1. Introduction 

The experimental realization of trapped quasi-one-dimensional atomic gases during the last few years 
[UlllinillllSllSllTllHlin] has provided a new impetus for the study of the effects of strong correlations 
on the physical properties of fundamental quantum mechanical systems of interacting particles. 
One very interesting aspect of this still emerging field, in contrast to traditional condensed matter 
physics, is that it is possible to follow the unitary time evolution in the absence of any dissipation 

[n mi mill mini US]. 

The absence of a general theoretical framework to investigate the non-equilibrium dynamics of 
extended quantum systems led to the study of many one-dimensional models whose dynamics can be 
solved exactly (see e.g. Refs. [ni[ia[Ii[20l[21J[2a[23[M[23[2i[2a[2i[2^ 
[35l [36] [37l [38l [39] 140] |4T] , but the list is far from being exhaustive) and inspired the developing of 
new numerical methods among which time-dependent density-matrix renormalization-group (DMRG) 
is the most successful |42| . giving the first opportunity to investigate the out of equilibrium behavior 
also of non-integrable systems as e.g. in Refs. I2Illl3lllliniliail7lllHllli[5ni[SI]- 

^From the study of specific models it is difficult to draw general conclusions on the non-equilibrium 
physics of quantum systems. Hence it is desirable to have frameworks that, at least in some relevant 
time windows, can give results valid for a large class of systems. In Refs. [24j [30] a method to tackle 
the time evolution after a quantum quench (i.e. a sudden change of a Hamiltonian parameter) has 
been proposed. Within this method the real-time path integral representation of a c?-dimensional 
quantum system is reduced to the thermodynamics oi & d + 1 dimensional classical system in a slab 
geometry. The time translation invariance is explicitely broken by the initial state that plays the role 
of a boundary state in the slab. When the Hamiltonian H governing the time evolution is at or close 
to a quantum critical point, one can use the renormalization group (RG) theory of boundary critical 
phenomena (see e.g. [S5]) to study the regime of time and length scales much larger than microscopic 
ones. In Refs. [24l|30| this approach has been mainly applied to the time evolution of one-dimensional 
(ID) systems since the 1+1 dimensional strip is described asymptotically by a boundary conformal 
field theory (CFT). 

Two very intriguing properties were derived by using this approach. First, connected correlations 
start forming only after two points are causally connected (horizon effect). Second, the large time 
asymptotic of correlation functions are the same as in a thermal state with effective temperature 
PcK = 4to (tq being related to the inverse of the mass gap in the initial state, see below) , despite the 
fact that the system is always in a pure state. These predictions have been confirmed by a number of 
exact calculations in specific models. This path-integral approach is generalizable to several different 
situations and it is very surprising, despite the success in the simplest case, that this has been done 
only in a few cases (in Ref. [33] for the so called local quench and in Ref. [53' for the evolution of the 
order-parameter in higher dimensional systems). 

In this paper we apply this imaginary time path-integral to the case in which the initial state 
is not translationally invariant but contains a domain-wall. We study the resulting time evolution 
governed by a gapless Hamiltonian using boundary CFT. A part from the per se interest, there are at 
least two additional reasons to study such an initial state. On one hand, the presence of a domain- wall 
generates a non-trivial transport, as shown by the exact solution in XX and XY chains [54]. On the 
other hand, in 2D classical critical systems, these are the typical boundary conditions used to set up 
the stochastic Loewner evolution (SLE), a mathematical rigorous approach to describe stochastic and 
geometric properties of 2D critical systems. Hence one must ask the question of whether SLE could 
be used to predict the time evolution of these quantum models. 

The paper is organized as follows. In Sec. [5] we review the imaginary-time path-integral approach 
to quantum quenches. We specifically show how to tackle the case of a domain-wall in the initial 
state by means of CFT. In Sec. [3] we apply the method to the Ising model for which all the relevant 
correlation functions in the upper half-plane have been calculated in Ref. [55]. In Sec. [4] we describe 
which features of this quenches can be extracted by general CFT scaling without the knowledge of the 
full correlators in the upper half-plane and then in the following section [5] these results are used to 
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give prediction for the time-evolution of the entanglement entropy. In Sec. [S]we argue about possible 
connections with SLE and we discuss the case of two domain-walls in the initial state. Finally in 
Sec. [7] we critically discuss our findings and open problems. In Appendix A| we report some useful 
formulas to perform the analytic continuations and in [Appendix B we derive the time evolution for an 
homogeneous quench of the two-point functions of different operators. 



2. The CFT approach to quantum quenches 

In Refs. [24l [30] it has been shown how to extract from path-integral the unitary time-evolution of a 
d-dimensional system prepared in a state IV'o) that is not an eigenstate of H . The expectation value 
of a local operator 0{r) at time t is 

{0{t, r)) = Z-i(V'o|e*''*-^^0(r)e-^^*"^^|i^o) , (1) 

where the damping factors e"*^^ have been introduced in order to make the path-integral representation 
of the expectation value absolutely convergent. Z — (?/'o|e~^'^'^|V'o) ensures that the expectation value 
of the identity is one. Eq. ([T]) may be represented by an analytically continued path integral in 
imaginary time, described by the evolution operator e~('^2~'^i)^, which takes the boundary values ipo 
on T — Ti = —e — it and on t — T2 ~ e ~ it. The operator O is inserted at r = 0. The width of the 
slab is 2e. ti and T2 must be considered as real numbers during all the calculation. Only at the end 
they are continued to their effective values ±e — it. The real-time non-equilibrium evolution of a d 
dimensional system is then reduced to the thermodynamics of a d -I- 1 field theory in a slab geometry 
with the initial state |'0o) playing the role of boundary condition at both the borders of the slab. More 
generally the operator O in Eq. ([1]) can depend on several points, and in the d — I case on which we 
focus below, we denote the set of these points, together with possible domain-wall positions in the 
initial state. 

We wish to study this setting in the limit where t and all separations \ri — rj \ are much larger than 
any microscopic length and time scales, so that renormalization group (RG) theory can be applied. If 
H is at or close to a quantum critical point, the bulk properties of the critical theory are described by 
a bulk RG fixed point (or some relevant perturbation thereof). Any translationally invariant boundary 
condition fiows to one of a number of possible boundary fixed points [SSI [57], and we may replace |'0o) 
by the appropriate RG-invariant boundary state |V'o) to which it fiows. The difference may be taken 
into account by assuming that the RG-invariant boundary conditions are not imposed ai t — ti and 
T2 but at r = Ti — To and t = T2 + tq. tq is given by the absolute value of the extrapolation length and 
it characterizes the RG distance of the actual boundary state from the RG-invariant one [52] . In the 
quantum non-equilibrium problem, tq is expected to be of the order of the correlation length in the 
initial state, that is the inverse energy gap [24l [30] . The effect of introducing tq is simply to replace e 
by e + Tq. The limit e ^ 0+ can now safely be taken, so the effective width of the slab results to be 
2to. For simplicity in the calculations, in the following we will consider the equivalent slab geometry 
between r = and r = 2ro with the operator O inserted at t = tq + it. 

If the initial condition (and so the boundary one in the imaginary time representation) is not 
translational invariant, the approach must be changed accordingly. For example, in Ref. 33J it 
has been shown how to deal with initial conditions corresponding to the joining of two semi-infinite 
ground-states of a critical Hamiltonian. The case of domain-walls in the initial states, in which we are 
interested in this paper, is also easily tackled with boundary CFT. In fact, as shown by Cardy [55] . 
any domain-wall at some position can be viewed as the insertion of an appropriate operator 
at position r^. The label a stands for the boundary condition to the right, and b to the left of the 
domain- wall respectively. $a6 is called boundary condition changing operator for obvious reasons. For 
any bulk CFT, there are a number of possible boundary states and consequently a number of boundary 
condition changing operators between them that are classified and their scaling dimensions are known 
(see e.g. |5^) in the most relevant cases. 

In the upper half-plane with a boundary condition having m domain-walls (including the point 
at oo) on the real axis, any n-point function can be obtained in terms of a (2n -I- TO)-point correlation 
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of the n operators plus m boundary condition changing ones. The 2n comes from the bulk points and 
their images. Increasing n and m technical difficulties to get these correlations increase considerably. 
For this reason, we focus on one- and two-point functions with two domain- walls (one at cxj). The 
difficulty of this calculation corresponds to that of four- and six-point functions in the complex plane. 

2.1. The general strategy by conformal mapping 

To obtain the time-dependence of a correlation function evolving from a state with one domain- wall (let 
say at r = for simplicity), we need to analytically continue the same correlation function calculated 
in a strip of width 2tq with boundary condition changing operators at r = on both the borders of 
the strip. The strip complex coordinate is w = r + ir with < t < 2to. The operators needed for the 
correlation functions are then inserted at t ^ tq + it. This geometry is depicted on the left of figure 
[TJ Note that we have set the sound velocity to unity. 



Z2= X2+ iy2 



2iTo 



wi = ri + iTi 



W2 = r2+ iT2 



- + r 




Figure 1. Left: Strip with domain-walls at = and w = i2To and points at r,; + ir, with t that 
must be analytically continued to tq -|- it. Right: Corresponding geometry on the upper half-plane 
after the mapping (|2]l. For the two-point correlation functions, the Euclidean distance between the 
two points is given by = (xi — x^)^ + (yi — ?/2)^- 



The strip geometry can be obtained by the conformal mapping of the upper half-plane z — 
with a domain-wall in the origin (and one at oo) 

7rt«/(2ro) 



1 



tanh 



4to 



■ly 



(2) 



We indicate in all the paper the upper half-plane with U and the strip with S. is the argument of z 
in U, as shown in figure [1] (right panel) . Note that the conformal map used here is different from the 
one used in Refs. [241 130j w{z) = (2To/7r)ln2; since we fix two marked points on the boundary (the 
two domain- walls) . 

In the case where O is a product of local primary scalar operators ^iiwi^Wi), the expectation 
value in the strip is related to the one in U by the standard transformation 



(3) 



where = 2hi is the bulk scaling dimension of <f>i. Often in the following we use the shorthand ^i(z) 
for $i(z,z). Note that, as usual, when comparing to lattice models, an eventual expectation value 
of in the ground state of Tl is supposed to have been subtracted off. The asymptotic real time 
dependence is obtained via the analytic continuation t — )■ tq -I- it, and taking the limit t, 3> tq. 
In the following sections we apply this method to several specific cases. 
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3. The Ising model 

We start our analysis from the Ising model for which several correlation functions (of the two primary 
operators spin and energy) in the upper half-plane with a domain-wall have been already calculated by 
Burkhardt and Xue [55| (using specific results for four- and six-point correlation functions previously 
derived [58j). In the beginning, we shall focus our interest to domain- walls created by changing fixed 
boundary conditions from a = -I- to 5 = — . We refer to this case as "H — initial conditions" . The 
initial position of the domain wall is at r = (unless noted otherwise) and a on its positive side. Then 
we separately discuss the case a — + and 6 = free, a domain- wall defined via changing from fixed to 
free boundary conditions. We will abbreviate this case by "-f/ initial conditions". All results here 
should be relevant for the time evolution of e.g. the Ising chain in a transverse field with Hamiltonian 
H = — \ Sii'^f *f+i + ^^i] (note that the state |-|-) here corresponds a spin ahgned in the x direction). 
This model is easily diagonalized through a Jordan- Wigner transformation followed by a Bogoliubov 
rotation, that map it onto a free fermion lattice model (see e.g. appendix of Ref. [20] )■ It is critical at 
h = 1 and hence described there by a CFT with central charge c— 1/2. 



3.1. H — initial condition: one-point functions 

Spin. The a expectation value with H — boundary condition in U is [55] 

, , „ . cos 9 , , , 

(4) 

where ^ is a non-universal amplitude, that can be fixed in a universal way through the ratio with 
the amplitude of the bulk two-point function [SSI ISO] (but we will never use the explicit value of 
A). This correlation function satisfies the correct (bulk) limiting behavior as y ^ -\-oo. Hence it is 
straightforward to read off that the scaling dimension of a is = 2ha = 1/8. 

The magnetization profile in the strip with the central domain-wall is obtained by using Eq. ^ 
and the conformal mapping Q 

/ NS ^ A,/o( JlX'"^ sinh(7rr/2To) 1 rr. 

^ \2tq) (cosh(7rr/To)-cos(7rT/To))i/2 [sin(^T/2ro)]i/8 • ' 

To extract the real-time evolution we analytically continue to t = tq + it, obtaining (r stands for the 
distance from the domain-wall) 

It w A fo ( ''^ \ ^^^ sinh(7rr/2To) 1 , , 

(a(r,t))=.AV y—j (cosh(7rr/ro) -hcosh(7ri/ro))i/2 [cosh(7rt/2To)]i/8 ' ^> 

which for asymptotic large time and distances tq <^ r,t reads 

(.,..0>=.4(-) ^^l^j^ (7) 

This result is depicted in the left panel of figure HJ 

Let us now discuss and interpret this first result. It is the consequence of two competing effects: 

• The states -I- and — are not eigenstates of the conformal Hamiltonian. Thus a half-line with -|- 
(or — ) initial condition evolves according to the results of Rcfs. [24l|3^. In particular, we expect 
{a{r,t)) (X e"'^*/^^'^", that is nicely confirmed by Eq. ([7]) in the case t < \r\. Far enough from the 
origin, the system does not feel the domain-wall in the initial state and the exponential decay of 
the one-point function is the same as for a translational invariant initial state. 

• The effect of the domain- wall (encoded in the first term of Eq. (jOJ)) arrives at the point r only 
at time t = \r\ (we recall that the speed of sound has been set to 1) and it is responsible for 
an additional exponential time decay e'^*^*"''"!-'/^'^''. This can be seen as a result of the quantum 
averaging of spins with opposite orientations. 

As usual [SU [30] , the effect of finite tq is only to smooth the time evolution close to the crossover point 
t^lrl. 
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Energy density. The time evolution of the energy-density operator can be obtained in analogous way. 
In the upper half-plane we have [55j 

,4cos2 6l-3 

(8) 



B- 



V 



i.e. the scaling dimension of e is Ag = = 1 (as well known for Ising) and there is a different 
non- universal constant B. Conformal mapping to the strip and the analytic continuation is worked 
out using the formulas in [Appendix A[ For t, r ^ tq , we obtain 

(9) 



{e{r,t)) 6 ( ^ ] e-'^*/^-" X 



t > |r|. 



See figure [2] (right) for illustration. The interpretation of the result is similar to the previous case, 
except that, contrary to the spin operator, the energy density of the -|- and — state are equal and so 
the result is symmetric with respect to the axis r = 0. For t < \r\, the time evolution is the same 
as in the absence of the domain-wall that only plays a role for t > \r\. However, for the energy, the 
domain-wall has a drastic influence, because the averaging of spins with opposite direction changes 
abruptly (at least for tq <C i, |r|) the value of the energy density. 

A cartoon of the time evolution of the one-point functions, e.g. in the transverse field Ising chain, 
goes as follows. The initial domain-wall state undergoes a spin-flip process which locally looks like 
I - +) ^ (I + +) - I - ->)/\/2. This process does not affect the two-site magnetization (in the x 
direction), but makes vanishing the single-site one, that before was maximum. Also the energetic 
balance is positive since the singlet state has lower energy than the product one in the critical Ising 
chain. Similar processes at later times t > I?"! are responsible for the negative local e and the suppressed 
magnetization (as in Eq. ([7]) inside the causal-cone of the domain-wall. We stress that this is only 
a cartoon of what happens locally because the real-time evolved state is highly entangled on a global 
scale (see the section about entanglement entropy). 
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Figure 2. Space-time dependence of one-point correlation functions from a H — initial condition. 
Left: (^(r, f)) for t = 20, 30 as function of r. Right (e{r, t)) for t = 18, 20 as function of r. Both the 
y axes are in arbitrary units and we fixed tq = 1. 



3.2. +- initial condition: two-point functions 

Spin-spin. Let us start with the spin-spin correlation function. In the upper half-plane U, we 
abbreviate a{zk) with Zk = Xk + iyk by cr(fc), fc = 1, 2. According to Ref. 55J, for boundary conditions 
H — , as shown in figure [1] the two-point function is given by 

((t(1)(t(2))" = —7- V"" + cos 01 cos H , sin 6*1 sin 02 ) 

V2 (2/12/2)^/** V Vu + u-^ J 

{ail))l_ {a{2))l_ ^ u^-l^ \ 

= ^yu-\-u MlH — ^ tan 01 tan 02 , (10) 

V2 V + 1 / 
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with the same ampHtude A as for the one-point function. Here 6^ are polar angles as indicated in 
figure [11 and u an auxiliary variable defined via 

„=(i + fM»)'". (U) 

Always using Eqs. ^ and Q we obtain the strip result S 

/^nW9^\S (^(l))+-(^(2))+- r-—T(. , "'-1 sin(7rri/2ro)sin(^T2/2To) ^ 
^ V2 V + 1 smh(7rri/2ro) sinh(7rr2/2To)/ 

where the Jacobians from conformal mapping have canceled out. The first factor does not contain any 
substantially new information since it is the product of two one-point functions. Analytic continuation 
of the various terms is performed in [Appendix A[ Putting everything together, we can write the time 
evolution of the two-point correlation function for i,ri, \ri — r2| 3> tq as 

({a{\)){cj{2)) 2t<\r,-r2\, 
H1H2)) = I (.(1)^.(2)) ^,(2.H..-^.l)/i6ro + ,^,^,.(2*-(|..|+k.|))/2.o) 2i > \n - r,\, ^'^^ 

with Si ~ signri, where m 1 in the first regime and w — > oo in the second (since 2t < \t\ — < 
I''! I + one can check from (|XTO)l that there is no third regime where small u — \ terms could mix 
with a large contribution from the s\Si term). This divides in different sub-regimes according to the 
relative values of ri, r2, \r\ — r2\/2, t. More precisely, we can distinguish two main regimes: 

• 2t < \ri ~ r2\'. the two points evolve independently since the horizon effect prevents from mutual 
interaction. The correlation decays as the product of the two one point functions and the connected 
correlation can be said to vanish (more precisely, its decay is subdominant). There are two possible 
sub-regimes: 

o / vr \ , f 1 t< min(|ri|, Iral) , 

WirM2»^A'{-) ...e-./-x|^„._„.„,„ ^^^J' ,14) 

the second regime exists only if |ri — > 2min(|ri|, |r2|) and we have written it in the case 
1^21 < (we will adopt this convention in this subsection). Let us denote la and lb these two 
subregimes. 

• 2t > \ri — r2\: One has then the following possible sub-regimes (which we can denote IIa,b and c, 
respectively) : 

- For |ri ~ r2\/2 < t < \ri\, \r2\ we find 

A2 / \ 1/4 

(a(l)a(2)) ^^(-] . (15) 



\/2 V^o, 

This is a quasi-equilibrium regime. Both points have space-like separation to the domain- 
wall, but time-like to each other, hence they equilibrate as if no domain-wall was present. 
Obviously this regime exists only if the two points are on the same side of the domain-wall 
S\S2 = 1. It only exists for a finite time interval, as eventually the domain- wall influence will 
be felt. 

For \r2\ < i < |ri| we have 

^ ^ ^ ^ ^/2\t^) I g.(t-|n|)/2.„ 2t> |ri| + |r2|, ^ ^ 

with an intermediate regime when correlation again evolves in time. Its absolute value first 
decreases down to the minimal value | (cr(l)cr(2)) Uin oc e-''(l''i-''2l+'*(l'^il-l''2l))/i6^o reached 
at time t = i\r\ \ + |r2|)/2. After that time it increases again. In addition, if the points are 
on opposite side of the domain- wall, there is a sudden change in sign (i.e. on scale tq) of the 
correlation at time t = (|ri| -I- |r2|)/2 — \r\ — r2\/2 hence only the second (increasing) part of 



Time evolution of ID gapless models 



8 



this subregime exists. Note that ahhough the absolute value vanishes on time scales smaller 
than To, it is possible to define a minimum value on scales much larger than tq, as can be 
seen on the two bottom figures in Fig. [3l 
— Very large times 2t > \ri — r2|, 2|ri|, 2|r2|, almost all terms cancel out and we may simplify 
to 

i.e. the correlation function is the same as in a thermal state with an effective temperature 
/?off — 4to as in the case without the domain- wall: the system has "relaxed" it. 




20 40 60 80 100 120 20 40 60 80 100 120 

t t 



Figure 3. Two-point function of tlie spin operator for four different quenches with (r2,ri) = 
(120, 90), (100, 30), (120, -90), (30, -60)). They are the exact results with tq = 1. When a single rj 
is negative (rj < 0) we plot the absolute value. Signs are indicated by the colors of the curves: black 
corresponds to + and red to — . 

Let us discuss these findings critically and describe which of these regimes are encountered as 
the system evolves. We start with the case of ri and r2 on the same side of the domain-wall, i.e. 
without loss of generality ri > r2 > 0. Then ri — r2 > can either be larger or smaller than 2r2. If 
(ri — r2)/2 < r2 then the evolution can be summarized as follow: 

• 2t < ri — r2'. the connected two-point function is constantly zero. 

• (ri — r2)/2 < t < r2: the two-point function has a plateau at e-^i^i-'^2)/i()To ^ Until this time 
the two points at ri,r2 did not yet feel the presence of the domain- wall. In fact the result is 
exactly the evolution found from a boundary state without domain- wall [241 130] , in analogy with 
the one-point functions. 

• At t = r2 , the part of the system under consideration finally enters in the horizon of the domain- 
wall and the two-point function starts evolving again with an exponential decreasing up to the 
time (ri -|- r2)/2. At this time, the correlation starts increasing toward the asymptotic value. 

• At t = ri the system reaches its asymptotic value and the two point correlation does not evolve 
anymore. The domain-wall has been relaxed. 
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An example for (ri — r'2)/2 < r2 is reported on the top left figure [31 where (ri,r2) — (120,90). The 
succession of regimes is thus la, Ila, lib, lie. 

The case 2r2 < ri — r2 turns out to be slightly different, and we give an example with 
(^i,''2) — (100,30) in the top right figure[3l In this case the sub-system feels the domain-wall before 
it is quasi-relaxed and so the succession of regimes is now la, lb, lib and IIc. 

When the two points are on different sides of the domain- wall, clearly | ri — r2 1 /2 (being the average 
of their absolute value) is always in between |ri| and |r2|, so that the plateau (quasi-equilibrium) regime 
Ila never exists. Hence the succession of regimes is la, lb, lib (second half only) and IIc. Let us point 
out an interesting difference to the previous cases: for short time the two-point correlator is negative, 
since siS2 = —1, and changes sign aX t = |ri — r2|/2 (when the systems goes from regime lb to lib) in 
order to reach a positive asymptotic value. This can be easily understood: in fact for 2t < \ri — 
each point is out of the horizon of the other and evolves keeping its initial sign. Connected correlations 
start forming only after t = \ri — r2\/2 when the two-point function changes sign. 




I I 



20 40 60 80 100 120 
t 




20 40 60 80 100 120 
t 



Figure 4. Absolute value of the energy-energy two-point function {t{r\,t)t{r2,t)) as a function of t 
for (ri,r2) = (120,90), and (100,30). We fix to = 1. Signs are indicated by the colors of the curves: 
black corresponds to + and red to — . 



As a general conclusion, although the correlation reaches the same asymptotic large time 
equilibrium value as if the domain- wall was absent, we found that it does exhibit an interesting non- 
monotonic time dependence, which shows that the domain-wall has a disordering effect on the system 
at intermediate times. A complementary picture of the various regimes is given below (see Fig. [11 



Energy- energy. We proceed to the two-point function of the energy-density operator e which has been 
calculated in Ref. |55| in the upper half-plane. We now map it to the strip and rewrite the resulting 
correlation function in a suitable way 

{e{l)e{2))\_^{e{l))\_{e{2))\_ 



1 



1) 



1-4 sin 



-1 1 -4sin2(6'i -^6*2) 



(l-4sin^ 6'i)(l -4sin^ 02) 



(l-4sin^ 6ii)(l-4sin^ 6*2) 



(18) 



Transport to the strip is simple and the analytical continuation of the various terms is presented in 
[Appendix A| Note that the factor — 1 indicates that the connected two-point function vanishes 
for 2t < \ri — r2\, as it can be seen from Eq. (jA.lOp . As above, it corresponds to the independent 
evolutions of the two points, which are not yet causally connected. Hence we focus on 2t > |ri — r2|, 
when the two points interact. In this regime m** — > 00 from Eq. (jA.lOp and one finds that only the 
second term in the parenthesis of Eq. (|18|) contributes. By analogy with the spin-spin correlation 
function, we find three different sub-regimes for 2t > |ri — r2|: 

(i) For t > l^il, |r2|, the system displays thermal like correlations and, in fact, we find 

2 



(£(l)e(2))+_=6M- 



^-7r|ri-r2|/2To 



(19) 



with the correct prefactor, and with /3cff = 4ro as it should be from general arguments (see below). 
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(ii) In the opposite limit t < \ri\, \r2\ (always with 2t > \ri — ^21), the correlation function is given by 
Eq. (jl9p . with the same prefactor. The two points are quasi-equilibrated since they do not yet 
feel the presence of the domain- wall. 

(iii) For \r2\ <t < \ri\ (we assume without loss of generality |ri| > \r2\) we find 

(e(l)£(2))+_ = -3B^ (^^^ ki-r.|/2ro ^ (20) 

where the result is the same as in the other two regimes, only being multiplied by —3. At variance 
from the spin correlation function, does not depend on time (apart from an obvious smoothing 
close to the changing times t — |r2|). This apparently strange phenomenon, visible in Fig. 
m can be traced back to the fact that {e{r,t)) for t > \r\ does not have any further exponential 
suppression, but only a discontinuity at |r|. 



Spin- energy. In the upper half-plane U, the spin-energy correlation function reads 



1 



1 



tan 01 sin 2^2 



(21) 



(.(l).(2))VH^(l)>V«2))Vx-^ 2 -4cos^.2-3^ 

The passage to the strip and the analytic continuation is easily done using formulae in [Appendix A[ 
Again for 2t < \ri — where u ~ 1, this two-point fimction is just the product of the one-point, in 
agreement with causality. 
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Figure 5. Absolute value of the spin-energy two-point functions (a"(ri , t)e{r2, t)) as a function of t 
for (r-i,r2) = (120,90), (90, 120), (100,30), and (-60,30). We fix tq = 1. Signs are indicated by the 
colors of the curves: black corresponds to + and red to — . 



Slightly more involved is the opposite case 2t > |ri — r2|, where asymptotically we have 

(a(l)6(2))+_ = (a(l))+_(e(2)>+_e'^(2*-l---l)/4- + ififHl^ . (22) 



Here Q{x) denotes the Heaviside function which is 1 for x > 0, and otherwise. We will not give a 
detailed discussion of the various sub-regimes. Instead we give some examples in figure [S] where we 
illustrate the time evolution for several choices of ri and r2. Amongst the various cases, the most 
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interesting and new feature when compared to (ccr) and (ee) is the very large time behavior: as t 
exceeds any other scale, i.e. for t> |ri|, |r2|, \ri — 7-21/2, we find that 

(CT(l)e(2)) + _ = 2yi6s2e-'^*(l/16+l/2)/rog-(|ri-r2|-2|r2|)/4ro^ (23) 

which decays to zero for t ^ oo. Indeed this feature is in agreement with the well-known fact that 
the finite-temperature equilibrium value of a two-point function of operators with different scaling 
dimensions vanishes (see also [Appendix B[ ). 



3. 3. +f initial condition: one-point functions 

Let us consider different boundary conditions: we suppose a -I- state on the half line from to +00. 
On the negative real axis the spins are left free to take either value -I- or — . 



Spin. In the upper half-plane the one-point correlation function of the spin is [S5] 

(a>^^-^2;-i/Vcos((?/2). 
Using results from [Appendix A{ the one-point function is asymptotically given by 



(a(r,i)) 



+/ 



A 



To 



1/8 



-7rt/16ro 



2-1/4 
e(r) 



\r\<t, 
\r\>t. 



(24) 
(25) 



See figure[6| (left panel) for an illustration. Note that the amplitude vanishes on the space like negative 
region r < —t as for homogeneous free boundary conditions, and in fact matches the result of Eq. ([7]) 
on the space-like positive side r > t. Inside the time-like cone it takes a non trivial value. 



Energy. In U, the one-point function for the energy is [55j 

e +/ = - cos 9, 
V 



(26) 



that (a part the constant and the scaling dimension) is the same as the spin expectation value with 
H — conditions and so we have the time evolution as 



+f 



To 



-TT(t-\r\)/2To 



\r\<t, 
\r\>t. 



(27) 



An illustration is given in figure [6| (right panel) . 
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Figure 6. Left; One-point function {cr(r, t)) with +/ initial conditions for t = 20 (blue) and t = 25 
(black). Right: One-point function (e(r, t)) for t = 20 (blue) and t = 21 (red). We fix tq = 1 and the 
j/-axis are measured in arbitrary units. 
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3.4^. +f initial condition: two-point functions 



Spin-spin. 
(a(l)a(2))+^ 



Again, 



Ref. 



[55 gives the half-plane result 
1 



V2 (2/i2/2)i/« 
(a(l))+^ (a(2))+^. 



v/(w + u-i) cos(6'i/2) cos(6'2/2) + (u - u-^) sin(6li/2) sin(6'2/2) 



1 



1 



tan(0i/2)tan(6'2/2). 



(28) 



V2 ' ■ V w^ + 1 

Relevant formulae for conformal mapping to the strip and analytic continuation are given in 
[Appendix A[ One finds similar regimes as for the H — initial conditions except that this two-point 
function is always positive, as a consequence of the non-negativity of one-point correlators. The main 
features are now standard. For time such that 2t < |ri — the two-point correlation function is just 
the product of the one-point ones, as m — > 1 in Eq. (|28p . For larger times, it is important to note that 
in the limit u — > 00 one should be careful in using the second line in Eq. ((28)) because it has the form 
X 00. It is simpler to use the first line properly mapped to the strip (as explained in [Appendix A[ ) 
that is always finite. For times such that 2t > \ri — \ one finds, using that u — *■ 00, a pseudo-thermal 
behavior with effective temperature PcS — 4to, for all three possible subregimes (noted IIa,b,c above), 
with 



(cr(ri,i)cr(r2,t))+/ 



1 



V2 



To 



1/4 



-7r|ri— r2[/16To 



(29) 



for both small and large time regimes t < |r'i|,|r2| and t > |ri|,|r2|. More surprisingly one finds 
the same time independent result (j29p for the intermediate time regime |r2| < t < |ri|, but with an 
amplitude multiplied by the factor 2~^/*. This is at variance with the behaviour of the H — case, and 
it is illustrated by two examples in figure [71 where two different pairs of values of ri, r2 are considered. 
We conclude that this kind of mixed domain-wall can be "dissipated" by the system, too. 

The other two-point correlation functions with (-1-/) boundary condition are known in the half- 
plane from Ref. |55j . We do not report their analysis here, since the main features of the Ising 
universality class are already present in the correlation reported so far. We rather highlight general 
features in the next section and show which information can be obtained from general scaling arguments 
only, without detailed knowledge of the full correlators. 
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t 



Figure 7. Spin-spin two-point functions (cr(ri , t)a{r2 ,t))^f as a function of t for (ri , r2) = (120, 90) 
(left) and (120, —30) (right). We fix tq = 1. When at least one r is negative, the two-point function 
is practically zero outside the two causal cones. 



4. General CFT results 

Some of the results previously obtained for the Ising CFT with particular boundary conditions have 
a direct generalization to an arbitrary theory with any domain-wall in the initial state between 
two boundary conditions a (to its right) and b (to its left). In fact any n-point correlation 
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(Ci ^ ^2)''^^"'' {'^i[zi) ■ ■ ■ ^n{znyfab '^^ primary operators satisfies the same diff'erential equation as 
the full plane 2n + 2-point correlation (zi) • • • $„(z„)$„(z„)<i>af,(Ci)*i'a&(C2)) and can be 

constructed from the same solutions, as discussed in full details in Ref. [55j . Ci and C2 stand for 
the start and end positions of the interface on the real axis, taken here to be (0,oo), and is the 
scaling dimension of the ab boundary changing operator $a6- Let us discuss general consequences of 
this property. 

The one-point function of any primary operator in the upper half-plane displays the scaling form 

($(z))l=A6y~^*/a*(cos0), (30) 

where A$ is the scaling dimension of <i> and all the dependence on the boundary state is in /*f,(cos0), 
with, as usual, cosO = x/{x'^ +2/^)^^^- In fact, /*f,(cos6') explicitly depends on the boundary CFT and 
is related to hypergeometric functions of the corresponding full plane 4-point correlation. 

A large class of primary operators displays a one-point function with a close analytic form (given 
in Ref. [55j). but in general they are quite complicated. However, from the previous examples, it is 
clear that we do not need the complete knowledge /^^(cosfl) to understand the quench. Indeed upon 
mapping onto the strip and analytical continuation we already found 

t<\r\, 

(31) 

Hence we need only the behavior for cos0 ~ 0, ±1. The case cos9 ^ ±1 means that in the upper 
half-plane the point z is far from the domain-wall and does not feel its effect. Consequently, f{cosO) 
takes the same value as for homogeneous boundary conditions a or b, related to cos 9 — 1 and — 1 
respectively, with associated amplitudes Aaa,bb = -4ab/*j,(±l). Conversely the case cos9 = regards 
any point inside the causal cone of the domain- wall (i.e. causally connected within the speed of sound). 
The amplitude /*^(0) can be either zero or finite, depending on the boundary CFT and the operator. 
We thus obtain 

when /* (0) ^ 0. Instead if /*(0) = 0, the decay becomes faster, as - f,-'^t/^^/^T,,-t,(^(t-\r\)/2Ta f^j. 
t > \r\ if /*^(cos0) ~ cos^ 9 around cos 6* — 0, as in the case of the spin correlation with H — boundary 
condition for the Ising model (with ^ — \). An even more detailed discussion of the one-point functions 
is reported in Sec. [6]in connection with SLE. In particular we show explictly in which cases 7*^(0) can 
vanish and that /*j(cos0) is analytic in zero, so in general = 1. 

The (relative) simplicity of these one-point function is related to the fact that bulk four-point 
functions can always be written as a general scaling part times a function of the unique cross-ratio of 
the four points. This is not the case for the two-point functions, which correspond to full plane six-point 
functions involving three independent cross ratios and the analysis is more complicated. However, we 
can give a general characterization. Considering a primary field $ such that ($(z)) 7^ 0, a scaling 
argument shows that we can always write the two-point function in U with boundary condition ab as 

(<f(zi)$(z2))l = ($(zi))^:,(<f(zi)):;,F*(u,cos^?i,cos02) , (33) 

where the function F^^ is (to the best of our knowledge) unknown in the general case, and related to 
full-plane six-point functions. The variable u is defined in Eq. pT|) is (one of the possible) cross ratio 
of Z1.2 and zi^2 (note that comparing with Ref. [30] one has — > 1 for u ^ 00 and 77 — > for u — > 1). 
Formula ([33]) still holds after mapping to the strip (replacing the index U by S) with the appropriate 
expressions for the variables u and cosfl^ repeatedly used above and in [Appendix A[ We thus find that 
the interesting regimes for the time evolution are encoded in particular limits of i^^. Namely: 



Time evolution of ID gapless models 



14 



• Short times t < |ri|, |r2| and points on the same side of the initial domain-waU. In this case, 
depending on the signs of r.i, we have cos 6i = ±1. This corresponds to two points extremely close 
to one boundary condition, i.e. 

i^*Kl,l) = F,,H, F^,{u,-1,-1) = FM, (34) 

where Faa{u) {Fbb{u)) corresponds to the two-point scaling function with homogeneous a (6) 
boundary condition. In this case the problem reduces to the one evolving from homogeneous 
initial condition. Note that A'^^Faa{u) = (u'*/4)^* for u — > cx3 as one recovers the bulk behavior 
and, upon mapping to the strip, thermal like behavior in the subcase i > |ri — r2|/2 (see also 
below). 

• Short times 2t < \ri — ri\, even on different sides. In this case u ~ 1, i.e. the distance of the two 
points is much larger than the distance from the boundary and -^^(1, cos ^i, cos ^2) — 1; meaning 
that the correlation function evolves like the product of the one-point ones. 

• Larger times 2t > \ri — r2\, the points are causally connected and u = e'^*^^*^''"!^''^!''/^'^'' '3> 1. 
Scaling requires that F^i^{u,rii,ri2) — m'*'^*9!'*{,(?7i, 772) in this large u limit, so that the time 
dependence factorizes in front and 

{Hr,,tMr2,t)),,^Al,(^^y^\-^^-\^^-^^\/^'^"^ (35) 

The very large time subregime t > |ri|, |r2| corresponds to rji = cos 9i = 0. Consequently, in U the 
two points are deep in the bulk and give the bulk two-point function, leading to the thermal-like 
behavior. The other sub-regimes, when t is in between the various "critical" times \ri\ can be 
extracted from ((35l) by considering the limits cos^i — > ±1,0 and may lead to additional time 
dependence. Note that AlJfMff,ir,2)<f>tbiVi,V2) ^ 4"^* as (771, 772) ^ (1,1) or (-1,-1). 

The same reasoning can be applied to any rt-point correlation function: for short times all the 
points evolve incoherently, whereas for very long times the correlation function is the same as in a 
thermal state at the effective temperature f3cS = 4ro. As for the homogeneous quench [30] it is easy 
to understand the technical reason that gives the effective temperature Atq both for asymptotic large 
times and also in the intermediate regime. As usual, finite temperature correlations can be calculated 
by studying the field theory on a cylinder of circumference P = 1/T. In CFT a cylinder is usually 
obtained by mapping the complex plane with the logarithmic transformation w — /3/(27r) log z, but 
one can equivalently use the inverse of Eq. (2) w = /3/7rarctanhz that differs from the logarithm 
for a mapping of the real axis. The form for the two-point function of a primary operator in the 
strip depends in general on the function F(u, cos 6i), but when we analytically continue, we find that 
u — > 00, i.e. the original points in U effectively are far from the boundary (more precisely their relative 
distance is much less than the distance from the real axis). Thus we get the same result as we would 
get if we conformally transformed from the full plane to a cylinder, and from the comparison of the 
two transformation the effective temperature is (3cS — 4ro. The same argument is easily worked out 
for the multi-point functions. 

5. The entanglement entropy 

Entanglement is a central concept in quantum information science. Moreover it is becoming a common 
tool to study and analyze extended quantum systems because of its ability in detecting the scaling 
behavior close to a quantum critical point ^61j . A powerful measure is the block entanglement entropy 
defined as follows. If the system is in a pure quantum state j^*), p — is the density matrix. 

Indicating with A a spatial subset of the system (such as a finite subset of spins in a spin chain) and 
with Pa — Ti'B P the reduced density matrix of the subset A, the entanglement entropy is just the 
corresponding Von Neumann entropy 

Sa = - TrA PA In pA , (36) 
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and analogously for Sb- When p corresponds to a pure quantum state, Sa gives a measure of the 
amount of quantum correlation between A and B, and one can prove that Sa = Sb- 
Sa can be calculated in a quantum field theory through the replica trick [62] 



(37) 

=1 

This is particularly useful in CFT because, when A consists of disjoint intervals with N boundary 
points with B we have that Tr transforms under a general conformal transformation as the N point 
function of a primary field <&„ with conformal dimension [62j Ajj — c/12(n— 1/n), where c denotes the 
central charge of the underlying CFT. The fields are usually called twist operators. In particular 
this implies that the entanglement entropy of a slit of length ^ in an infinite system is given by [63l [62] 



Trpl = c„ - ^ iog--cl, (38) 

\a / 6 a 

where a is an UV cutofi^ (e.g. in a spin chain is the lattice spacing) and c'^ a non universal constant. 

Clearly within this mapping to correlation function we can calculate the time-dependence of the 
entanglement entropy. For the case of homogenous quenches this has been done in Ref. [50] , finding 

— t<i/2 , 

SA{t) ^ { (39) 



t > 1/2 



12to 

This result can be easily explained in terms of the causal scenario longly discussed there, and it is 
due to pairs of coherent quasi-particles emitted from any point and reaching one the subsystem A 
and the other B. In the case of A and B corresponding to two semi-infinite lines is a straightforward 
generalization of the previous result. In fact, we have only to consider the one-point correlation 
function of a twist operator ($„) that in the upper half-plane is ($„(2:,z))u = c„[Imz]~'^" (where c„ 
is the analogous of c„ for the boundary CFT introduced in [S2]). Mapping to the strip, analytically 
continuing, and taking the limit of large time, directly gives that the entanglement entropy grows 
indefinitely with time as 5a (i) = 7rct/12To. Let us mention that this indefinite grow of entanglement 
is the main reason why numerical simulations (e.g. a la density matrix renormalization group) do not 
perform efficiently for this kind of time evolution [64j . 

The natural question arising is how this result is modified by the presence of the domain-wall 
at r = 0. Let us start with the case of two semi-infinite lines with A = [r, oo) and B = R\A its 
complement. In this case the entanglement is given by the one-point function ($„(r, i)) of a twist 
operator at position r. From the general result reported above Eq. ([5^ . we know that 5a (i) increases 
linearly with time for t < \r\, as expected because the time evolution is the same as in the absence of 
the domain-wall. For t > \r\ we also find from Eq. (j32p a linear increase of the entropy, but its rate in 
principle could change if (0) is vanishing. However, the properties Ai = 1 and Tr p = 1 trivially 
leads to AabfabW ~ ^- Analyticity in n and cos 6* near zero makes then /*f,"(cos6' 0) non zero 
close to n = 1 and this ensures that the domain- wall does not change the linear increasing rate of the 
entanglement entropy, but only changes it up to an additive constant. The finite change, with respect 
to e.g. a homogeneous initial conditions of type a, is given by Sab = -c^n(-^afc(/a6" (0) ~ fabi^))\n=i, 
independent of r with |r| < t. Hence the entanglement of the two parts being initially separated by 
the domain- wall increases linearly with only a universal additive constant Sab distinguishing from the 
homogeneous quench. This constant is computed below, in the particular cases relevant for SLE. 

For the entanglement entropy of a finite slit A — [ri,r2] with £ = \ri — r2\, general features are 
simple consequences of the scaling form in Eq. p3|) . As long as t < |ri|,|r2|, the time evolution 
remains the same as in the absence of the domain- wall (i.e. initial linear increase and then saturation). 
Conversely for large times t > l^ij, |r2|, £/2 the entanglement entropy saturates to a value proportional 
to £. In the intermediate regimes, that depend on the nature of the domain-wall, from the general 
scaling again only finite universal shifts compared to the homogeneous quench can be possible. 



Time evolution of ID gapless models 



16 



6. Connection with SLE 

Geometrical and statistical properties of a two-dimensional critical system in the half-plane with an 
interface starting at the boundary are nowadays widely studied by means of the so called stochastic 
Loewner evolution (SLE) that allows for a rigorous mathematical derivation of some CFT results (see 
e.g. Refs. [5S] as reviews). We do not enter here in the beautiful technical details of the approach, but 
we will limit to state some results and ask whether they may have some connection with the quench 
problem considered here. An SLE interface receives a CFT interpretation as being created by inserting 
a boundary changing operator $12 at the boundary. In some cases, such as Ising, it corresponds to 
a domain wall and an ah boundary change. Hence for simplicity we stick here to the Ising notations, 
keeping in mind that in this section we deal with boundary condition changes relevant for SLE^ at 
generic value of the parameter k. These may be quite different from Ising, and are detailed in e.g. 

Refs. [eg. 

One of the most celebrated result concerns the rigorous derivation of Schramm's formula |66| . 
Starting from the boundary the interface explores the full upper half-plane and ends at infinity. 
This curve (called the SLE trace) divides the half plane in two parts. Schramm's formula gives the 
probability that one point lies to the right of this curve. In the continuum limit, this probability only 
depends on the angle 9 that the point makes with the boundary (that is the same 9 we used in all the 
paper). Schramm's formula gives this probability as 



satisfying the "boundary conditions" p{9 = 7r/2) = 1/2 (reflection symmetry) and p{9 = tt) —0 (point 
at the negative real axis) . The parameter k is the diffusion constant of the underlying stochastic process 
and turns out to define different universality classes (e.g. the Ising model corresponds to k = 3) and 
is related to the central charge. 

Since SLE focuses on geometrical properties of interfaces, one could try to push the continuation 
to real time i, and define a space-time interface geometrically using the path integral formulation. This 
can be done in principle by appropriate continuum limit (in the time direction) of lattice models such 
that the interface remains well defined. Consider the Ising case as a typical example. Specifying an 
interface means specifying the local spins -I- (to the left) and — (to the right) along a simple curve in 
the r, i plane (since it wanders this entails creation and annihilation of domain- walls as t increases). 
One must then superpose the complex amplitudes associated with all configurations satisfying that 
dynamical constraint. These are continuations of the real positive probabilities used to construct SLE. 
Schramm's dimension zero operator additionally constraints the interface to pass right of a space time 
point (i.e. to be connected to one of the two boundary conditions a or 6 without crossing the interface). 
Clearly this qualitative picture needs to be made more concrete. 

Can p{9) be used to determine some properties of the quench problem studied here? In order 
to address this question, let us map p{9) to the strip according to Eq. ([2]), and then perform the 
analytical continuation r — > tq -I- it. A simple calculation gives 



that is a positive number smaller than 1. Hence it could in principle still be interpreted as a probability. 
In the limit t, \r\ 3> tq, using the properties of the hypergeometric function, we find 



independently from k. This result is suggestive and could be thought, in a sense which remains to be 
made precise, as a probability of being in the state evolved from +, that is clearly zero for r < —t and 
1 for r > t. A plateau develops (see figure [5]) for r < \t\ and by mirror symmetry it can here only take 
the value 1/2. It would be interesting to understand how this plateau forms from the path integral 
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Figure 8. Probability p(r, t) as given in Eq. I|41|l for different times t. Here, the parameters are 
tq = 0.1 and k = 3 (Ising model). At large times, a plateau around r = at p = 1/2 emerges. 



formulation, presumably from interference effect. From the point of view of a real time relativistic 
theory it is however rather natural. 

Are there other probabilities to which this could be compared? A priori the quantum state 
evolution is fully deterministic and apparently there is no probability into the game. A possible way to 
introduce a probability proceeds as follows. Consider, in a lattice model, the reduced density matrix 
Prit) of the degrees of freedom localized in r (i.e. for a spin chain the spin in r). Clearly for r > t, 
Pr{t) equals the analogous one obtained from the evolution of a homogenous + boundary condition 
Prit)- Analogously for r < — i, Pr{t) is p^it)- An interesting possibility is that the reduced density 
matrix in the center region, — r < i < r, is simply a combination of the above two, in the present case 
by symmetry it can only be Pr{t) — {Pr (t) ~^ Pti^))/'^- interesting question then is whether the 
1/2 found from Schramm's formula is the same factor appearing in front of the p^{t) reduced density 
matrix. It is easy to check in the Ising case, from our result for the averaged spin, that this could only 
hold strictly for infinite time (the leading decaying term in e^'^^/^'^o does not agree). Unfortunately, 
in that limit, the factor 1/2 being imposed by symmetry does not provide a test. Tests would require 
more complicated situations, e.g. translating the results of Refs. [S71 [S5l [59] in terms of quantum 
evolution. For instance, in Ref. [67j the probability Pmiddie{cot 9) was obtained, that a point z in the 
upper half plane lies in between (left of one, right of the other) two interfaces starting at nearby points 
on the real axis, conditioned to the fact they do not cross (i.e. going to infinity). This may be relevant 
to determine whether, in presence of two domain-walls in the initial state, e.g. an interval [0, vq] of — 
spins in an otherwise homogeneous -I- state, the point r ^ ro at time t still feels the presence of the — 
in the initial condition (however, this provides only part of the answer, since one should add the event 
where the domain-wall annihilates at late times, but this also can be obtained, in principle, within 
SLE). The large time continuation yields Pmiddie (cot 6) D^, the constant defined in Ref. [67], for 
all t > \r\, which has a non trivial value, e.g. £'k=3 = 0.921866 for Ising. 

To close on a less speculative note, let us give the explicit form for the one-point function of any 
primary field of scaling dimension A in presence of the boundary conditions which are relevant for 
SLE. It is known that the boundary changing operator which creates an SLE„ interface is $12, with 
conformal dimension A12 = 2/112 — (6 — k)/k. The result of Ref. [55] is then Eq. (|32p with 

/afc(cos 9) = Ci 2F1 (a, p, 1/2, cos^ 9) + C2 cos 9 2F1 {a + 1/2, /3 + 1/2, 3/2, cos^ 9) , (43) 

where for convenience we have absorbed the amplitude Aab into fat, with parameters 

a, = i (1 - 2Ai2 ± [(1 - 2Ai2)2 + 12A(1 + A12)] '^') = («; - 4 ± ^ {k - 4)^ + 8kA))/{2k) . (44) 
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We determine the constants Ci and C2 from the boundary conditions /af,(±l) aX — tt (respectively 
to the right /left of the domain- wall) , that yields 

^ _ 1 Iab{l) + /ab(-l) ^ _ 1 /afe(l)-/afc(-l) . 

' 2 2i^i(a,/3,l/2,l) ' ' 2 2Fi(a + l/2,/3+l/2,3/2,l)- ^ ^ 

Using the identity 2^1(0,6,0,1) = T{c)T{c - a - h) / {T{c - a)T{c - b)), we find c-a~h= (8- k)/2k. 
Hence both denominators are finite for < k < 8 (the denominator associated to Ci equals 1 for all k 
for A = and this is consistent with the result for the twist operators with Ai = in the entanglement 
entropy section). However we find that Ci has a pole for A = A^/^^ — f + 5) ("^ + « ~ 5) ^'^'^ '-^z 
for a1^^ = ^{n+ 1) (n + |) both for non-negative integer n (note that the latter vanishes for A = 2 
for all k). Hence for such operators the symmetry must be imposed for the result to be meaningful. 
Obviously one recovers the Ising results Eqs. @ and ^ for k = 3 inserting respectively A = 1/8 
(spin) and A = 1 (energy). 

Using the continuation formula (|3ip one finds formally the large time asymptotics for t > \r\ is 
given by Eq. with 

AabfabiO) ^ Ci + C2e-'^(*-l'-|)/2^» + . . . , (46) 

where the second decaying term becomes the leading one when Ci = 0, else it is a rapidly decaying 
subleading correction. 

Note that the Schramm formula is recovered, for some choice of boundary condition, considering 
an operator with scaling scaling dimension A = 0. One gets then (3 = and a = 1 — A/k and 

fabicosO) = Ci + C2 cos 61 2^1(1/2,3/2 -4/k, 3/2, cos^ 61) = Ci + C2 cot 612^1(1/2, 4/k, 3/2, - cot^ 6*) , 

using properties of hypergeometric functions. It remains to be seen whether a dimension zero primary 
operator, relevant for the quench, can be identified. In that case it should be related to Schramm 
formula (or be trivial). 

We can now compute, in the case of an SLE interface, the universal additive contribution Sab to 
the entanglement entropy as discussed in Section [5l Inserting A = c/12(n— 1/n) which yields the 
appropriate values for a„ and /3„, and using the symmetry /*f,"(^l) = .f*fc"(l); ^'^'i 

Sab = -dM{an,(3nA/2,l)]-\=i = c "^^^ ^o^^J^^^""^ ' (47) 

3(k — 4) 

where ipix) stands for the Polygamma function and, in terms of k, the central charge is c = 
(6 — k){3k — 8)/(2k). Notice that for k < 4 (and in particular for the Ising model) this universal 
shift is negative. 



6.1. Two domain-walls 

Here we consider the case of two domain- walls in the initial state. It is more general, but in the context 
of SLE, it corresponds to a two interface problem and can be studied, for appropriate boundary 
conditions, by insertion of two $12 operators on each side in the strip. The corresponding strip 
geometry is depicted in Fig. [9l However, the complete correlation functions mapped on the upper 
half-plane need the insertion of four changing operators, that are in principle calculable but very 
cumbersome. To make progress, we first consider the limit of a small b domain in a line a when results 
can be obtained from the fusion of two <I'i2 operators, which produces $13 operators. At the end, we 
finally argue that for asymptotic large time the result is valid for any finite fixed size domain b, not 
necessarily small. According to physical intuition this means that no matter how big is the b domain, 
it is always small compared to the infinite remaining a part. 

Let us consider the operator O = Yii^ii'^i^^i) being a product of primary operators. When 
l^ci — X2\ = \x3 — x^l = ^ is small, taking into account the insertion of boundary condition changing 



Time evolution of ID gapless models 



19 







aba 



w = r + ir 



^\\\\\\\\/ / / / / JI^WWW^ 

a b a r 

H 1 h 



~e/2 +i/2 

Figure 9. Strip geometry with domain-walls at ±£/2 and ±£/2 + 2jro. 



operators, and the (generic) fusion algebra $i2<&i2 ~ 1 + ^13, we obtain 



4 



\ k=i / s 

= {0)l + {{O $13(0))^. + {O <^i3{2iTo))l) + ^2Ai3 (0 ci>i3(0)cI>i3(2zTo))^ + . . . , (48) 

where the dots denote the contribution of descendant operators, which are easily seen to include 
analytic 0{P) corrections as compared to the leading term which corresponds to homogeneous 
boundaries. The dimension of the $13 operator is, in the SLE context: 

Ai3 = ^ . (49) 



In terms of SLE of each term in Eq. p5|) can be interpreted as follows. The second term represents 
the contribution of two interfaces starting both very close to zero (respectively to 2iTo) on the same 
side of the strip and which do not immediately annihilate near their origin (but rather annihilate in 
the bulk). The last term corresponds to the possibility where the two interfaces cross the whole strip. 
It has thus some general relation to the calculation of Ref. [67] , relation which can be made precise for 
a dimension zero operator. We recall here that we use Ising notations for convenience only but that 
we consider generic k and the appropriate boundary changes |65j). 

Let us now specialize to a single primary field O — ^(w) of conformal dimension A, being inserted 
at w = r + ir. As explained above, at the leading order in these can be obtained simply from ($)^, 
($$13(0))^, and ($$i3(0)<i>i3(2iro))^. To calculate this one should map the result in the upper half- 
plane, that are obtained again by fusing two $12. The calculations are straightforward. For example, 
applying a Moebius map to Eq. (|30p one obtains the one point average in presence of a single finite 
domain, ($(z))^,„ = ^''^^=($(^)$i2(0)$i2(0)u = 2/-^/*fc (cosV') with cos^- = Rei^ z{l - z)/z{l - z)) 
(the constant Aab is absorbed in /*j again). This gives the lowest order terms in Eq. P5|) . with the 
correct amplitudes, from the expansion of /*;,(?7 = 1 — C)//*b(l) = 1 + CaC^^^^^ + BaC + ■ ■ ■ for small 
( = ^{e/y)^sm'^e + 0{e'^). it yields in particular ($(z)$i3(0))u ^ Caismef'^^^y-^-^'^ . Putting 
everything together we finally obtain 

{<^)ltaba - (<f>"a (^1 + Ca (^)^" (siu 0)^^13 + ^ (^^^ + {^)'^'^ 5*(cOS0) + ..^ , (50) 

with Ba - a/3/(l/2 + a + (3) ^ 4A/(8 - 3k) and Ca = 2^^! (1/2 - a, 1/2 - /3, 1/2; 1) {Ci/fab{l)) + 
2F1 {1 - a,l - P, 3/2; 1) (C2//ab(l)), where Ci and C2 are given in (gSl), a, (3 in ([44]), and the relation 
given below (|45|) to ratio of Gamma function should be used to define the relevant limits. The general 
expression for the function g'^{cos9) is given in Ref. [55j in terms of hypergeometric functions, using 
the fact that $13 being degenerate at level three, ('I>(z)<l'i3(0)<f>i3(oo))^ obeys a simple differential 
equation. Note that the descendants will give non trivial corrections to ([SU)) . but these are always 
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subdominant to the leading correction and, quite often subdominant even to subleading ones. Hence 
in the general expansion in small {i/y) at fixed cos0, the leading correction comes from the {£/y)^^^ 
term only for A13 < 2, i.e. for SLE type boundary conditions, k > 8/3. For k < 8/3 the Ba term 
dominates. Note however that even for k > 8/3 the Ba term dominates if the amplitude vanishes 
Ca = 0. This is the case for the Ising model k = 3 if <i> is either the bulk spin or energy operator, since 
fabiv) is analytic in these cases as discussed in Section S) The last term is often subleading, although 
in principle there may be cases, for k > 4, where Ca vanishes and the leading correction is given by 
the last term. Another example consists in the 0{n) loop model when each interface carries a different 
0{n) index hence they cannot annihilate. The general case depend on the operator $, on the model, 
and on n, and can be read from the above results. 

Mapping to the strip and continuing to the real time quantum dynamics, we know that for t > \r\ 
the half plane variable {£/y) sin^ 9 in ([50|l maps to ^e^'^^*^''-'/^'^'', a variable which is small both in the 
small £ and the large time limit. Hence we find for large times i > I?'! 

(^(r, t))aba = ($(r, t))ail + C„ (^e-^f*-'')/^-" + i?a£2e-(t-r)/r„ _^ _ ) _ (5^) 

where one of the two correction terms dominates depending on whether A13 is larger or smaller than 2, 
and of whether Ca vanishes or not, as discussed above. The scaling variable governing the correction is 
then ie~'"*^^'^°, so no matter how large £ is, provided it is finite, for long times the leading correction is 
always of the above form. Physically this corresponds to the property that asymptotically the b domain 
is washed out and the time evolution is the same as from an homogeneous a initial condition with only 
exponentially small corrections. Although we have not attempted any calculation it is physically clear 
that any initial state with a local domain structure immersed in an infinite a initial state will behave in 
a similar fashion, possibly involving a spectrum of higher Ai^2Ar+i exponents, corresponding to 2A'^ SLE 
curves. Applying to A ^ A„ and following the arguments of Section O one also easily sees that the 
difference in entanglement entropy as compared to the homogeneous case decays exponentially in time, 
with the same decay rates as in ([^ and prefactors given by —dn{Ca,Ba)\„=i- Exphcit calculation 
shows that the leading amplitude does not vanish. Similarly, in presence of an odd number of domain 
walls in a finite region in the initial state, the difference in entanglement entropy as compared to the 
homogeneous state, is expected to converge at large time, with similar exponentials, to the value of 
Sab given by (gT]). 

7. Conclusions 

We presented a detailed study of the time evolution of a ID quantum system evolving according to 
a gapless Hamiltonian from an initial state with one domain-wall. All our findings can be easily 
interpreted in terms of quasi-particles excitations emitted at time t = 0, and then propagating trough 
the system at a finite speed given by the sound velocity Vs — 1- According to this scenario any 
correlation function has different time regimes. For example, in figurefTUlthe time evolution of the two- 
point correlation functions has different regimes resulting from the competition of the quasi-particles 
emitted from the middle point and those emitted from the domain-wall. The same interpretation is 
also valid for the entanglement entropy. 

The importance of this simple physical picture is that it easily allows to understand limit and 
generality of the CFT results. In fact, the behavior outside the causal cones is only consequence of 
the finite velocity of sound and so it must be generally true also for gapped systems, as pointed out 
in details in Refs. [24l[30]. Oppositely, lattice models, even in the thermodynamic limit, have a full 
spectrum for the velocity of excitations, due to non-linear dispersion relations. These slow quasi- 
particle give corrections inside the causal cone to the results just derived (see Ref. [30] for a careful 
discussion of this issue). 

As in the case of an homogeneous quench, we showed that for large times all the correlation 
functions are the same as in a thermal state with effective temperature f3eS = 4to. This means that 
the presence of the domain-wall does not affect the asymptotic state. In fact, this feature does not 
come unexpected since the energy of this state is only slightly larger than the translational invariant 
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one. Following this interpretation, we conclude that a state with an infinite number of domain-walls 
should change the effective temperature. 

We finally discussed the possibility that the geometric picture obtained from SLE may be continued 
from imaginary to real time. Whether this may lead to new information for quantum quenches is not 
yet clear at this stage, but has been critically discussed in the present paper. We have also obtained, for 
the type of domain wall considered in SLE, the universal change in entanglement entropy as compared 
to an homogeneous initial state. Finally we considered the case of few domain walls in a finite region, 
and obtained the leading time decay of one point functions. 




4-+++++++++++++++++++++++++++++++++ 



Figure 10. Causal interpretation of the results for all the two-point functions (0(ri, i)0(r2, t)). In 
the case shown with (r2 — ri)/2 < ri < r2 and > 0, there are four different time regimes, shown 
with different colors in the figure. 1) For 2t < r2 — ri, the two points evolve independently (blue 
region). 2) For (r2 — ri)/2<t<ri, the two-point function is quasi stationary because there is not 
yet an influence of the domain- wall (green region). 3) For r\ < t < r2, the domain- wall enters in the 
dynamics (violet). 4) For t > r2, stationary behavior with quasi-thermal correlations (red). 
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Appendix A. Useful formulas for the analytic continuation 

The strip geometry is obtained by the conformal mapping of the upper half-plane z = x + iy to the 
strip w = r -\- ir given by Eq. It is useful to write it explicitly in cartesian and polar coordinate. 
We have 

g™/(2ro) f-KW 

= tanh 



sinh(7rr/2To) 



Rez = 
Im z — 



cos(7rr/2To) -I- cosh(7rr/2To) ' 

sin(7rr/2To) 
cos(7rT/2To) -)- cosh(7rr/2ro) 



2 ^ cosh(7rr/ro) - cos(7rT/ro) 

2[cos(7rT/2To) +cosh(7rr/2To)]2 " ^ 

Upon continuation t ^ tq + H one has sin(7rT/2ro) — > cosh(7rt/(2To)) and cos(7rT/ro) — > — cosh(7ri/To), 
but cos(7rT/(2To)) continues to an imaginary value. Hence the denominator must disappear in any 
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physical quantity, and it does, being compensated by the Jacobian 

' ^ ^' 4to' ^ ^' 2ro(cos(^r/2To) + cosh(^r/2To)) ^sin(7rr/2To) ' ^ 
The polar angle 6 is mapped to the strip according to 

„ X \/2sinh(7rr/2To) . „ y v^sin(7rT/2ro) 

cos = — = — ■ — - — i — '— — -, — — — ■ , sm — — 



\z\ [cosh(7rr/To) — cos(7rT/ro)]^/^ ' \z\ [cosh(7rr/To) — cos(7rr/To)]^/^ ' 

After analytic continuation we have 

^/2sinh(7rr/2To) . f 1 * < \r\, 

cost) ^ — ; ^ — ■ -"TTT- ~ Sign r X < ,^ ,,.,„ 

[C0sh(7rr/To) +COsh(7rt/To)]l/2 ^ | g-7r(t-|rl)/2ro i>|r|, 

sin^ ^ v^cosh(^t/2ro) ^ f e-(*-M)/2ro ^ < |,|^ 

[cosh(7rr/To) + cosh(7ri/To)]i/2 |i t > \r\. 

tan^ ^ 2^2^111^ ^ sign, e-it-H)/2ro ^ (A.4) 
smh7rr/2To 

where here and in the following ~ (and in the main text) stands for f, r ^ tq. 

We will also need the following functions of the angles that are trivially obtained from the previous 
formulae 

sin26l =2sin6icos6i~2signre-''l*-l''ll/2^°, (A.5) 
e /l + cos^V''' Z^^"") 

I TT fl t>\r\, 

^ e /I -cos 61 , , 

tan- =\ ~< ef-r) (A.7) 

2 V l + cos^ 1 lli-lZ t<\r\. ^ ' 

K B(r) 

4cos^^-3= ( 8sinhVr/2ro -1 + 2 sign (|r| - ^) (A.8) 

\cosh(7rr/ro) — cos(7rT/ro) / 

Here, G(r) denotes the Hcaviside function which is 1 for r > 0, and otherwise. Note the divergence 
for r < — t in tan^/2. However, this divergence cancels with the one-point function with free boundary 
condition and we do not need the first order correction in tq that makes it finite. 

For the two-point functions, we need to study the combination of different variables. In particular 

we need 

%Ly2 \ ^ {xi - xif + (j/i + yif ^ cosh(7r(ri - r2)/2ro) - cos(7r(ri + T2)/2to) 
(? ) {x\ - x-if ^ {y\ - y-if cosh(7r(ri - r2)/2To) - cos(7r(Ti - r2)/2ro) ' 
that is easily analytical continued to T1/2 — > tq -|- it as 

^4 _ ^ ^ cosh^(7rt/2ro) ^ [ 1 + e"(^*-l''^-''^l)/^"° ^ 1 In - ra] > 2i » tq, 

" COSh2(7r(ri - r2)/4To) \ eT(2t-|ri-r2|)/2ro ^1 2i > In - rsl > TO. 

In particular all the terms like (u^™ ^ 1) vanish for 2t < \ri — r2\ and the ratios (u^™ — Ci)/{v?"^ + C2) 
tend to 1 for 2t > |n — r2|, for all constants Cj. 
For combinations of polar angles, we have 

5^e^"(*-l''il'/2^'' ±S2e-^(*-l''^l)/2^° kiMrsI < t, 

sin(6'i ±6(2) = sin6iicos6»2±sin6i2Cos6'i ~ { S2 ± si e-''(l''=l-l''il)/^^° |n| < t < \r2\, 

s-^ e^it-l'-^D/^^o ± S2 e'^(*-l''2l)/2^o |n|, |r2| > t, 
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siS2e-''(*-l'^il)/2^°e-'^(*-|'-^l)/2^'' T 1 \r1\Ar2\ <t, 
008(6*1 ±612) = cos 6I1 cos 02 T sin 6*2 sin 6*1 ~ { siS2e~'^(*~l''il)/2^° T e''^*"'''''^/^^° ki| < t < \r2\ 

S1S2 T e"(*-l'-i"/2rog-(t-k2|)/2ro i^^i^ i^^i > 

where again Si — signs^. 

Appendix B. The two-point function of different operators in the quench from a 
homogeneous initial state 

In Refs. pil [5U] it has been shown that when the initial state is a translation invariant state, the 
long-time correlation functions are the same as a pseudo-thermal state with /3off = 4ro. Then it follows 
that the correlation of operators with different scaling dimensions must vanish, but we still do not 
know how they fall off to zero. In this appendix we fill this small gap. 

We first consider the energy-spin correlation in the Ising model starting from a homogeneous 
initial condition (we only consider fixed, since free boundary conditions gives a correlation identically 
zero). In the upper half-plane U the correlation function is [70] 

-I- 1 

{a{l)e{2))l = {a{l))l{e{2))l^^ , (B.l) 

where (ct(1))!J_ oc y^^^*, and (e(2))'|_ oc y2^, and u is defined as in pT|) (see [Appendix A] for an 
explanation how to manipulate it). For 2t < \ri ~ r2\, we have u 1 and the correlation function is 
the product of the one-point functions, in agreement with the causal interpretation. Oppositely when 
2t > |ri - r2|, because ^ g7r(2t-|ri-r2|)/2To^ ^lave 

(CT(l)e(2))+ CX e-'^*/16^0e~''*/^^"e^(2*~l''i"''^l)/*^° = g-^t/16rog-7rlri-r2|/4ro ^ (-g_2) 

i.e. the two-point function is going to zero with an exponential rate. 

Another simple example is the Gaussian model defined by the continuum Hamiltonian H oc 
J^^^CV (t>{z))'^ dxdy . The field (j) is not primary, but its exponential e*""^^^-* with arbitrary a is primary 
with scaling dimension c? 12. The general correlation function is given by 

^gm</.(2i)g-i/3,/,(5:2)^ ^ ga^G(2;i,2:2)-"Q''G(2i,2i)/2-/32G(z2,Z2)/2 (-g 3-) 

with 0(21,22) = ((/'(zi)0(z2)). Because of the restriction to the upper half-plane we have the Green 
function 

G(zi, Z2) = Goo(zi - Z2) - Goo{.zx - 2.2) (B.4) 

where Grx.[z) — 6*00(0) — log \ z\\s the Green function in the plane (we adsorbed the standard tt factors 
in the definition of 0) and we assumed a fixed boundary condition (free ones are pathological because 
the one-point function is identically zero). We then have 

/ ia</<(zi) -i/3</<(z2)\ _ \ -/rinl"'' nr u^"^ fB ^"l 

(e e )- |^^_|„2/2|^^_|/3V2 1^12/^12! «^„2/2^^^V2" • 

Using the previous formulas for the mapping to the strip and the analytic continuation we have for 
2t < \ri — r2\ the independent evolution of the two points, while for 2t > \ri — r2| 

^gia0(zi,t)g-*/30(z2,t)^ _ 

g-7rQ^t/4rog-7r/5^t/4rog7r(2t-|ri-r2|)Q/3/4To _ g-7r(a^+/3^ -2a/3)t/4To g-7rQ/3|ri -rj | /4ro j-g 

that always decays exponentially to zero because (a^ + (3'^ — 2a/3) > 0. 

It is easy to show that this indeed is a general feature of the homogeneous quench for different 
operators for any initial condition. In fact the, two-point function of different operators $1 and $2 in 
the upper half-plane can be always written as (when ($1) ^ identically) 

($i(l)<i>2(2))" = (<i>i(l))"($2(2))"i^(«^) , (B.7) 
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where F{u'^) can be related to the hypergeometric function determining the four-point functions in 
the bulk [71j. However its complete expression is not required: we only need the asymptotic behavior 
of ~ 1) and F{u'^ ^ 1). From general scaling arguments, we have F{1) = 1 (i.e. close to the 

boundary the two-point function factorizes in one-points ones) and F{u'^ S> 1) ~ u^^^*" . The exponent 
Xb can be calculated for any pair of operators and any boundary condition from the bulk four-point 
functions [71]. Here, we only stress that if $1 = $2 then Xb — x\~ X2 to recover the right behavior of 
the bulk two-point function, and in general Xb < {xi +X2)/2 for unitarity. Putting everything together 
we have 

g-^a;it/2rog-^X2t/2ro 2t < Iri - I , 

^-T!Xit/2To ^-TTX2t/2To ^irXi(2t-\ri-r2\)/2Ta 2t > \ri — r2\ 

Note that when the operators $1 and <i>2 have the same scaling dimensions the time dependence 
disappear and we get the usual thermal like correlator. Oppositely, since xi + X2 — 2a;f, > 0, the 
correlation decays to zero as anticipated. 



($i(ri,t)$2(r2,0)> 
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